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Abstract 

Ultra high energy (UHE) photons can initiate electromagnetic showers in magnetic field. We 
analyze the two processes that determine the development of the shower, e + e~ pair creation and 
synchrotron radiation, and derive formulae for the angular distribution of the produced particles. 
These formulae are necessary to study the three-dimensional development of the shower. 
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I. INTRODUCTION 



Magnetic fields play a foundamental role not only for the acceleration and propagation of 
charged cosmic rays, but also for the absorption of neutral particles, photons and neutrinos, 
if the magnetic field is sufficiently strong or the particles have sufficiently high energy [|TJ. 
In particular, photons can initiate an electromagnetic shower in magnetic field which is 
analogous to the showers produced in matter; the main features of such showers (longitudi- 
nal development and particles energy spectra) has been analyzed in Refs. [^|, [| under the 
assumption that momenta be orthogonal to the magnetic field. 

Important examples of electromagnetic showers in strong magnetic fields are the 7 and 
radio emission in pulsars (3| and blasars || in active galatic nuclei 0. Strong magnetic 
fields of the order of 10 12 G can be found in the proximity of pulsars 0, || |, (Tcfl : in 
magnetic field of this order of magnitude photons loose energy by radiating photons (photon 
splitting |T|) or creating e + e~ pairs ]12], [13| which feed the cascade by producing more 
bremsstrahlung photons. In such fields even neutrinos radiate photons |14], [15|, |I6|1 , create 
e + e _ pairs |L6], [17], or even W particles |JI9| . 

If the energy of the primary particle is sufficiently high, an electromagnetic shower can 
develop even in weak fields, such as those present in the interstellar medium or in the 
vicinity of stars and planets |2(], |21[ . A very important such a case is the shower produced 
by UHE (E > 10 18 eV) photons in the magnetic field outside the earth atmosphere: this 
early shower influences the later atmospheric shower. These UHE photons are predicted 
by top-down theories as possible explaination of the experimental spectrum of ultra high 
energy cosmic rays |22], |23], ^4], |(| ) . 

In this paper we study the two processes that are the building block of the electromagnetic 
shower in magnetic field: synchrotron radiation by UHE electrons and e + e~ pair production 
by UHE photons. In particular, we derive formulae for the angular distribution. The precise 
dependence of the angular distribution from the field strength and from the energies of 
the particles is necessary to determine important features of the phenomena, such as the 
three-dimensional development and the lateral spreading. 

A foundamental question is whether the three-dimensional development of the shower 
could experimentally discriminate between UHE air showers originated from a primary pho- 



ton or from a primary proton (or heavier hadron) [27], 28], |29| and, therefore, discrim- 
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inate between competing theories of the high energy tail of the cosmic ray spectrum. In 
fact, UHE photons start the shower well outside the atmosphere producing an additional 
lateral spread to the subsequent atmospheric shower relative to a shower originated by a 
proton. In addition a shower that begins outside the atmosphere is less affected by the 
Landau-Pomeranchuk-Migdal effect 0, [H], 0, 

Another context where it is important the precise knowledge of the angular distribution of 
the particles produced in the electromagnetic shower is the modeling of the pulsar emission. 
For instance in the polar cup model proposed by Sturrock H, [7|, || high energy electrons, 
due to the intense (~ 10 12 G) magnetic fields, follow the field lines to minimize synchrotron 
radiation energy losses: photons are emitted in a narrow cone nearly parallel to the field 
lines. The contribution of pair production to the photon interaction length depends on 
the magnetic field component orthogonal to the photon momentum: since momenta of the 
particles in the shower are nearly parallel to the field, the precise emission angle might be 
critical for the shower development. 

In the following Section |H] we introduce the notation, derive formulae for the synchrotron 
radiation (magnetic bremsstrahlung) by UHE electrons/positrons, discuss these formulae 
and show some relevant plots. In Section |TTTj we make the analogous derivations and discus- 
sion for e + e~ pair creation by UHE photon in magnetic field. The last Section |IV| is reserved 
to our conclusions. 



II. SYNCHROTRON RADIATION 



Synchrotron radiation (or magnetic bremsstrahlung) from ultrarelativistic electrons has 



been studied by many authors, see for instance Refs. |L3], [34], |35], |36|, |37fl , where many 
results mediated over the angular distribution can be found. For our study of the angular 
dependence of both for the synchrotron radiation and the pair production we shall follow the 
approach of Berestetskii-Lifshitz-Pitaevskii-Landau (BLPL) ||35|| . In the following discussion 
we shall assume that the electron momentum p is perpendicular to the magnetic field H: 
eventually we discuss results for the general case in the last section (conclusions). 
We recall some of the relevant notation. The characteristic parameter is 

- H E m 
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where H is the stationary magnetic field, m and E the electron mass and energy and H c is 
the critical field 
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m c 



Hc = = 4.414 • 10 13 G , (2) 

en 



which is a natural quantum mechanical measure of the magnetic field strength ||34|| . The 
four-momenta of the incoming and outgoing electron and the one of the emitted photon are, 
respectively, P = (E/c, p), P' = (E'/c,p f ) and K = (hu/c,k). In most of the formulae we 
shall use h = c = 1. 

We are interested in the limit when both the ingoing and outgoing electron are ultrarel- 
ativistic, E 3> mc 2 , and the field is relatively low H -C H c . In this limit the quantization 
of the electron motion is not important, since E and AE = Tiuo are much larger than 
Too® = hceH/ E, and we can use the semiclassical approximation of BLPL with the electron 
following its classical orbit. 

The differential probability to emit a photon is [)35| : 



with 



/+oo 
dtV fi (t) 
-oo 



2 d 3 k 

(2vr) 3 



V fi (t) = e I d 3 r r f & ^ A* = -^—{f \Q{t)\i)e wt 



were ipij are the wave functions of the initial/final electron state, A is the vector potential, 
and a are Dirac matrices; Q(t) is the corresponding time-dependent Heisenberg operator. 

Using the completeness relation and introducing the variables r = £2 — 1\ and t — (p2 + 
ti)/2, the probability of photon emission per unit time is: 

dw e 2 d 3 k /,+0 ° 



dt uj (2ti) 2 



dr(i\Q + (t + T/2)Q(t-T/2)\i)e- iWT , (3) 



where d 3 k = u 2 dujdVL = u 2 duj sin ddddip with -d the angle between k and p, and (p the angle 
between H and the projection of k on the plane orthogonal to p. 

In Eq. (Q) the probability can be expanded in powers of r, since the main contribution 
comes from small values of r, when there is superposition between the amplitudes. In fact 
the values of r for which the superposition is significant can be evaluated using semiclassical 
arguments. For kinematical reasons ultra-relativistic electrons radiate in a narrow cone 
"0 ~ m/E: the amplitudes for the emission sum coherently along a section of the classical 
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electron path where the direction changes of an angle $ ~ m/E, i.e., uqt ~ m/E, giving a 
formation time 

1 H r h H r 



m 



m 



1.29 • 10" 21 (^ c /fl')sec" 1 . 



(4) 



Elu eH m H mc 2 H 
To leading order in ujqt ~ m/E the resulting non-polarized photon emission probability is: 



dw 



e 2 uo r+°o 



dtdujdVt 



(27T) 



x exp 



dr 



E 2 + E 



f2 



, 2 mr 
4£' 2 - M + EE' 



(5) 



-i — out 1 — n ■ v , 
E' V 24 



where v is the initial electron velocity. 
The integration in dr of Eq. (pi) yields 



dw 



where 



dtduodQ, n 



(E 

X = 2 — 
\m 



£ 2 + £' 
EE' 



i2 



-X - 



2/3' 



i(X) 



2£ 2 / \kE' 



k is the characteristic parameter defined in Eq. (|l|) and $(Jf) is the Airy function |38 



1 f°° ft 3 

®(x) = - dt cos xt H 

Jo V 3 



(6a) 



(6b) 



(7) 



In the ultra- relativistic limit the differential probability of photon emission, Eqs. (E|), 
depends only on the angle d between the electron and emitted photon directions, threfore 
dVL = 2vrd(cos$). 

It is often useful, for instance when writing the cascade equations, to express the result 
in terms of the fractional energy carried by the photon u = uo/E; then equations (@) become 



with 



dw 2a 
dt dftio d cos $ h 



X -I- 

m 



l + (l-u) 
(l-«) 



-X 



u 



m 
2E 2 



cosfi 



'1 — u)k 



u 



2/3" 



d>(X) 



E 2 



1 + 2(1 - costf) — - (1 - COST?) 



1 — U)K 

u 



2/3 



2/3 



(8a) 

(8b) 
(8c) 



m* j \ (1 — u)k j 

where a is the fine-structure constant. In the ultrarelativistic limit it is convenient to use 

J2 _ , I 77A2 



2(1 — cos??) ~ $ ~ (m/E) as angular variable. 



Equations or, equivalently, Eqs. @ are our main new result, together with the anal- 
ogous result for pair creation, Eqs. (|TTD, result which is presented in the next Section, i.e., 
the angular dependence of the produced particles. 

The integral of Eq. (|8|) in dcos$ yields the differential emission probability per unit 
energy: 

where 

u \2/3 / u \ 2/3 



^w) =Ur^yi (9b) 



The result in Eqs. (0) agrees with previous calculations 13! 

Other important quantities in the study of the electromagnetic shower in magnetic field 
are the differential energy loss for unit time and unit of photon energy: 

dE a, 9N9 / f°° , , , l + (l-«) 2 $'(0\ , s 

---{mc 2 ) 2 u[ §{x)dx + —± - ^ , (10) 



dtdu h \J£ ' 1 — it £ 

where £ is given in Eq. (9_b), and the total energy loss per unit time (emissivity): 

dE a, 9 ^k 2 f°° 4 + 5kx 3/2 + 4k 2 x 3 . , 

— = -— mcT — / ; ^— x$'(x)dx . (11) 

Note that, as expected from the lack of other dimensional scales in the limit of E ^> m, 
the spectral emissivity in Eq. ([TOD depends only from k and u (scaling), and the emissivity 
in Eq. ( O ) depends only from the characteristic parameter k, apart from the overall factor 
(mc 2 ) 2 . 

Before studying the angular dependence of the emitted photons let us discuss the total 
energy loss as a function of k, Eq. (|TT|) , and the energy dependence of the differential energy 
loss, Eq. (03). For k < 1 the energy loss goes to zero as k , therefore we limit our discussion 
to the more physically important 

The main feature of Eqs. (§) and (|T0D which determines both the energy and the angular 
dependence is the presence of the Airy function <$>(x) that goes to zero exponentially for 
large values of x. Only for the purpose of this discussion we use x < 1 as a simple threshold 
value (note that $(l)/$(0) = 0.381, while $(2)/$(0) = 0.0984), i.e., we assume for the 
purpose of this discussion that most of the photons be emitted for value of x ^ 1 (the 
discussion does not change if we use 2 or 3 intestead of 1), and use the ultrarelativistic limit 
2(1 - costf) ~ d 2 ~ (m/E) 2 < 1. 
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According to this criterion the differential energy loss per unit of photon energy, Eq. (fLOD , 
is large when £ < 1, and, therefore, most of the photons are emitted with a fractional energy 
u that verify the condition 

"*TT- K • (12) 

i.e., photons with a large fraction of the electron energy are emitted only for relatively 
large values of k. In addition, since the energy loss is proportional to the energy fraction u 
carried out by the photon, the energy loss goes to zero with u. Figure [I] summarizes this 
discussion by showing the probability of emission as function of the energy fraction with an 
arbitrary normalization of 1 at u — for different values of k. Curve 1, which corresponds 
to the largest value of k = 100, has a peak for values of u close to 1, while curve 6, which 
corresponds to the smallest value of k = 0.5, is peaked at values of u close to 0.1 and goes 
to zero at u = 0. 

The same criterion applied to the angular distribution, i.e. the constrain 
implies that the angle within which most of the photons are emitted depends on the energy: 



1 = VT/l^l (14) 




u 



when £ < 1. If instead £ > 1 the angular distribution decays exponential from the value 
dEjm = with a width proportional to a/ k(1 — u)/u. In other words photons with large 
energy fraction, u ~ 1, can be emitted within a smaller angle than photons of low energy 
u< 1. Figure (|2|) demonstrate this effect by showing the probability of emission as function 
of the variable v/2(l — cos $)E/m ~ dE/m for four choices of the couple {u, k). Dashed 
courves (1 and 3) have k — 0.1 and show how the angle becomes more than an order of 
magnitude larger going from u = 0.5 (1) to u — 10~ 3 (3). The same effect is show for k = 1 
going from curve 2 to curve 4. The same angle grows more than a factor of two going from 
k = 0.1 1 (3) to k = 1 2 (4). 

III. PAIR PRODUCTION 

The amplitude for pair production can be obtained from the amplitude for synchrotron 
radiation using the cross-channel symmetry p5[. Results mediated over the angular distri- 
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bution can also be found in |U| [3C], |34], |3^, |39], [40], [41 . 

The calculation follows closely the steps in the previous Section with the necessary formal 
differences. The characteristic parameter \ is analogous to k in Eq. (fl|) with the substitution 
of the incoming-electron energy E with the incoming-photon energy hu: 

H hu 



X 



H r mc 2 



(15) 



we use a different notation for clarity. Now the four-momenta of the incoming photon and 
of the outcoming electron and positron are K — (lu, k), P_ = (EL, p_), and P + = (E + , p + ). 
As in the previous section, we work in the limit when incoming and outgoing particles are 
relativistic and the field is relatively low, H <C H c . 

After performing the cross-channel transformations (E,p) — > (— E + , — p + ), (E',p') — > 
(E_, p_), (uj, k) — ► (—uj, — k), and substituting the emitted-photon phase space with the one 
of the created positron (or electron), we obtain the analogous of Eq. fl5|) for the probability 
of pair production by an unpolarized high-energy photon that propagates othogonal to 
the magnetic field (the probability is summed over the final spins and integrated over the 
electron direction, if we measure the positron energy: the roles of electron and positron can 
be exchanged). 

This formula can again be expanded in the formation time t/, Eq. (fj), and the result is: 



dw 



dtdE+dQ„ 



e 

7T 



El 



LU 2 E 



\XuE_ 



2/3' 



where 



A =21 — 1 I I ! I 



m 



cos'd 



$(A) 



2/3 



16a) 



(16b) 



my V V 2EU~~~"J \ X (E+E-) 
Again there is no dependence on the azimuthal angle (p and dQ + = 2nd cos d. Introducing 

the fractional energy carried out by the positron, v = E + /u, Equations flUf ) become 

2/3" 



dw 



2a 



dtdE+d cos "Q 



v 2 + (1 - v)' 

a-v) 



■X + v 



v(l - v)x 



$(A) 



(17a) 



with 



X 



E 4 



m 



m 



2Ei 



cos$ 



Ei 



1 + 2(1 - costf)^- - (1 - cosi?) 



m 1 - 



v(l - v)x 
1 



2/3 



2/3 



(17b) 
(17c) 
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and 2(1 — cos??) ~ d 2 ~ (m/E) 2 the convenient angular variable in the ultrarelativistic limit. 

As in the case of the synchrotron radiation, we can cross-check our result by integrating 
Eq. ( |TT| ) over the polar angle and comparing the differential emission probability per unit 
energy with Refs. 0, |3§ [J§ ^ |T| : 

*» -..W-fr^.^) (18 a) 



dtdE + \uJ \J 6 y ' E + E_ £ 

where 

w 2 \2/3 / ! X 2/3 



(18b) 



X E + E -J \xv(l-i 
The differential emission probability per unit time and unit of fractional energy v is: 

im _.w ^ r tWttc _^+(i-»)'^ (19) 



cftdt! 7i Tio; yj^ t>(l — v) £ 

with £ defined in Eq. ( |18b| ), while the total pair production probability is: 

dw a{mc 2 ) 2 f°° ( 2x*l 2 + l/x 



dt h hu 



£,,. { J?J. + X)v> ) w ■ (20) 



Note again that the differential emission probability depends only from x and v and the 
total pair production probability only from x (scaling), apart the dimensional constant in 
front (mc 2 ) 2 , as a consequence of the lack of other dimensional scales in the limit of E ^> m. 
Note also that the pair production probability, differently from the synchrotron emission, is 
exponentially suppressed in the limit of x ^ 4, since the Airy function decays exponentially 
for large values of its argument. The physical cause of this suppression is the presence of a 
threshold for pair creation. In the following discussion we consider the range x ~ 4. 

As in the case of the synchrotron radiation the main feature of Eqs. ( O) and (|T9| ) that 
determines both the energy and the angular dependence is the exponential suppression of 
the Airy function $(x) with growing x; we use the threshold value x ^ 1 and work in the 
ultrarelativistic limit i?~m/£<l. In addition the pair creation probability is symmetric 
in the two variable E + and E_, i.e., is symmetric respect to the point v = 1/2. 

In this case the chosen criterion imples that the differential energy loss per unit of elec- 
tron/positron energy, Eq. flT9[) , is large when £ ;$ 1, and, therefore, most of the e + e~ are 
emitted with a fractional energy v that verify the condition 



1 1 ^ 1 ^ A 1 
4 X 2 V4 x 
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when x > 4, i.e., e + e~ can be emitted with one of them carrying a large fraction of the 
photon energy only for relatively large values of \. Figure [5] summarizes this discussion 
by showing the probability of emission as function of the energy fraction with the area 
arbitrary normalized to 1 for different values of \. The smallest value of \ = 0.1 (curve 
1) is strongly peaked at v = 1/2, while the largest value of x — 100 (curve 6) has a much 
flatter distribution with peaks at values close to v = and v = 1, in spite of the fact that 
the distribution must go to zero at exactly v = and v — 1. 

The angular distribution is described by exactly the same constrain found for the syn- 
chrotron radiation 

•^( 1 + (^) 2 )(^) 2/3 - 1 ■ (22) 

which implies that the energy-dependent angle within which most of the pairs are emitted 
is: 

^ < yV(i - v) x f 3 - 1 = y/i/F^ (23) 

where we are considering only values of £ < 1. If instead £ > 1 the angular distribution 
decays exponential from the value §E/m = with a width proportional to a/(1 — v)vx- In 
other words pairs that share equally their energy v ~ 1/2 can be emitted within a larger 
angle than pairs where one of the two particles carries a large part of the energy, v ~ 1 or 
v ~ 0. Figure (f|) demonstrates this effect by showing the probability of emission as function 
of the variable 2(1 — cos n 9){E /m) 2 ~ $ 2 ; the distribution has been arbitrary normalized such 
that it is equal to 1 at dE /m = 10 _1 . Going from an asymmetric distribution of the energy, 
v = 0.1, (curves 1 and 3) to a symmetric one, v = 0.5, (curves 2 and 4) the angle becomes 
wider; the same happens when x increases from 4 (curves 1 and 2) to 100 (curves 3 and 4). 

IV. CONCLUSIONS 

In this paper we have studied the angular dependence of photons emitted by UHE elec- 
trons and of the e + e~ pairs emitted by UHE photons in a static magnetic field: this de- 
pendence is needed for detailed studies of the electromagnetic cascade in magnetic fields, 
such as those initiated by UHE cosmic rays in the geomagnetic field, or by charged particles 
emitted by pulsar. 
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The main results are shown in Eqs. (§) for the magnetic bremsstrahlung and in Eqs. (|17D 
for the pair production. For simplicity we have sketched the derivation of these formulae in 
the case of propagation in the plane orthogonal to the magnetic field, but it is possible to 
show, as we have explicitely verified by performing the appropriate Lorentz transformations, 
that the same formulae are valid in the general case if H is substituted with H±, the 
component of the magnetic field perpendicular to the propagation. 

These results are also plotted as function of the angle for different values of the fractional 
energy in Figs. |^ and |] and briefly discussed in text. The angle becomes wider for larger 
values of the characteristic parameter and for smaller energy fraction (synchrotron radiation) 
or more symmetric energy fraction (pair production). 

We have verified that our results, when integrated over the emission angle, reproduce the 
known results for the differential in the energy and total probability of emission. 
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FIG. 1: Differential synchrotron-radiation emission probability as function of the fractional energy 
carried by the photon u = lo/E, Eqs. dlO|), for several values of the characteristic parameter 
k = (H/ H c )(E/m): curve 1 (2, 3, 4, 5, 6) corresponds to k = 100 (10, 5, 2, 1, 0.5). All curves are 
normalized such that the area is 1. 
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FIG. 2: Angular distribution of the synchrotron-radiation emission probability, Eqs. (||), for two 
values of the fractional energy carried by the photon u = uj / E and of the characteristic parameter 
k: curve 1 corresponds to u = 0.5 and k = 0.1; curve 2 to u = 0.5 and k = 1; curve 3 to u = 10 -3 
and k = 0.1; curve 4 to u = 10 -3 and n = 1. The angular variable is ^2(1 — cos"d)E /m ~ $E/m; 
the normalization is such that curve is 1 at i? = 0. 
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FIG. 3: Differential pair-production emission probability as function of the fractional energy 
carried by the electron/positron v, Eqs. (|l9|), for several values of the characteristic parameter 
X = (H/H c )(w/m): curve 1 (2, 3, 4, 5, 6) corresponds to % = 0-1 (0.5, 1, 4, 10, 100). All curves 
are normalized such that the area is 1. 
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FIG. 4: Angular distribution of probability of emitting a positron (electron) with fractional energy 
v by a UHE photon, Eqs. (fl7|) , for two values of v and of the characteristic parameter X = 1: curve 
1 corresponds to v = 0.1 and x = 4; curve 2 to v = 0.5 and x = 4; curve 3 to v = 0.1 and x = 100; 
curve 4 to v = 0.5 and x = 100. The angular variable is \/2(l — cosi!})E + /m ~ i9E + /m; the 
normalization is such that curve is 1 at i? = 0. 
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